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Abstract.—Global warming and the decline in precipitation threaten wetlands worldwide, and lakes in some 
regions are in the process of drying. Amphibians, since they are water-dependent, will be the creatures 
most affected by the rapid habitat losses due to climate change. Especially for amphibian species which are 
endemic, the situation will be more serious in terms of its impact on biodiversity. Therefore, in this study, we 
determined the climate characteristics specific to the habitats of an endemic amphibian species, Pelophylax 
caralitanus. According to the Representative Concentration Pathways (RCP) climate change scenarios of the 
ICPP, we analyzed whether the climatic characteristics specific to these habitats will exist in 2050 and 2070 
under the criteria of RCP 2.6, RCP 4.5, RCP 6.0, and RCP 8.5. The results are quite alarming for Pelophylax 
caralitanus. According to the RCP climate change scenarios, the climatic conditions in the present habitats of 
this endemic amphibian species will not remain stable in that the potential habitats in Southwestern Anatolia 
will be dramatically reduced and the appropriate habitats of P caralitanus around the Turkish Lake District will 
completely disappear, while some new potential habitats will emerge in the Northwest Aegean region of Turkey. 
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Introduction 


The Beyşehir Frog, Pelophylax caralitanus, is an 
endemic species of Turkey which is distributed across the 
Mediterranean region of Turkey called the Turkish Lake 
District (Arikan et al. 1994; Ayaz et al. 2006; Baskale et 
al. 2017; Dusen et al. 2004; Kaya et al. 2002). According 
to the literature on the distribution of P. caralitanus, it 
inhabits permanent wetlands with rich aquatic vegetation, 
including permanent ponds, rain pools, streams, rivers, 
and irrigation channels. Pelophylax caralitanus prefers 
abundant vegetation around lakes and ponds as habitats 
(Başkale and Çapar 2016). In particular, emergent 
vegetation present in the water and the presence of 
bushes and weeds around the water body constitute an 
ideal area for the release of eggs and seeking shelter from 
predators. Pelophylax caralitanus 1s listed by IUCN as 
Near Threatened (NT) because of ongoing threats from 
habitat loss and overexploitation (Oz et al. 2009). Since 
it has a relatively wide distribution and presumed large 
populations, it is unlikely to experience declines that are 
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rapid enough to qualify it for listing in a higher Red List 
category. 

Climate change is one of the greatest problems 
threatening biodiversity and ecosystems in the 21 
century, leading to the extinction of many species 
(Sinervo et al. 2010; Walther et al. 2002). The habitat 
preferences of amphibians are closely related to 
temperature, precipitation, and wetlands, so amphibians 
are one of the most sensitive animal groups to changes in 
climatic conditions (Enriquez-Urzelai et al. 2019; Ortiz- 
Yusty et al. 2013). Temperature was found to be positively 
correlated with the detection probability of P. caralitanus 
(Başkale and Çapar 2016). Wetland ecosystems, ın which 
amphibians complete their life cycles, are among the 
areas most affected by climate change, which creates 
inhospitable conditions for amphibian reproduction 
(Desta et al. 2012; Hopkins 2007). According to the 
literature, the losses of amphibian habitat will be 
high due to the increase in annual temperatures in the 
Mediterranean basin, and it is predicted that there will be 
new gains in the north (Araujo et al. 2006). In addition 
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to the difficulties amphibians face from possible climatic 
changes, the pathogens that cause amphibian deaths, 
such as Batrachochytrium dendrobatidis in the southern 
part of the world, have been detected with increasing 
frequency in the north due to climate change (Berger et 
al. 1998; Cohen et al. 2019; Erismis et al. 2014). 

The IPCC (Intergovernmental Panel on Climate 
Change), which has been investigating climate change 
for 30 years, presented new climate change scenarios 
under the name of “Representative Concentration 
Pathways” (RCP) in its Fifth Assessment Report (IPCC 
2014). These scenarios estimate that global warming 
will change atmospheric CO, concentrations by 2100 to 
levels of 1,370 ppm for the RCP 8.5 scenario, 850 ppm 
for RCP 6.0, 650 ppm for RCP 4.5, and 490 ppm for 
RCP 2.6. These scenarios predict that temperatures in the 
world may increase by anywhere from 1.5 °C to 5.8 °C 
(IPCC 2014). 

It is important to know the climatic requirements of a 
species in order to predict how that species will respond 
to climate change. The MaxEnt approach and software 
have been widely used in recent years to determine the 
current climatic demands of a species and to predict 
changes in that species’ habitats in the future (Hendrick 
and McGarvey 2019; Kıraç and Mert 2019; Untalan et al. 
2019; Zhang et al. 2018). MaxEnt gives more accurate 
results with less data in smaller areas compared to 
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Fig. 1. Turkish Lake District and species presence data. 


Amphib. Reptile Conserv. 


35°0'0"E 


THs 


other methods, e.g., DOMAIN, BIOCLIM, and GARP 
(Hernandez et al. 2006; Wisz et al. 2008). In addition, 
MaxEnt allows categorized and continuous data to be 
processed together (Phillips et al. 2006; Phillips and 
Dudik 2008) and it creates habitat suitability maps in 
addition to the outcome outputs (Elith et al. 2006, 2011; 
Hernandez et al. 2006, 2008). 

The aim of this study is to investigate the ecological 
niche of P caralitanus and evaluate the consistency 
and variations in the predicted potential distributions 
of this endemic species in the Turkish Lake District in 
the Southwest of Turkey under the IPCC’s four future 
climate scenarios (RCP 2.6, RCP 4.5, RCP 6.0, and RCP 
8.5). The results will help us to understand how this 
endemic amphibian species will experience the effects of 
climate change in relation to distributional shifts. 


Materials and Methods 


Study site. The Turkish Lake District (Fig. 1) is located 
in Southwest Turkey and includes the western part of the 
Taurus Mountains. This region contains many different 
tectonic water bodies such as Beyşehir Lake, Eğirdir 
Lake, Burdur Lake, Akşehir Lake, Salda Lake, and 
Acıgöl Lake. The Turkish Lake District is a transitional 
region between the Mediterranean and continental 
climates (Baskale et al. 2017; Dusen et al. 2004). 
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Data collection. Presence data for P caralitanus were 
recorded from 89 sites in the Turkish Lake District area 
based on our field studies performed between 2010 
and 2016, as well as previously published papers (Le., 
Ayaz et al. 2006; Baskale and Capar 2016; Baskale et 
al. 2017; Dusen et al. 2004; Kaya et al. 2002) (Fig. 1). 
The sampling sites and their coordinates are given in 
Appendix 1. Sampling from habitats during this study 
was approved by the General Directorate of Nature 
Conservation and Natural Parks in the Ministry of 
Agriculture and Forestry, Turkey (protocol number 
18.10.2010/61288) and the Animal Care and Use 
Committee at Pamukkale University (protocol number 
20.09.2010-PAUHDEK-2010/021). Individuals of the 
target species were identified in their natural habitats 
and were not exposed to stress. No measurements or 
experiments were performed on the species. 


Bioclimatic data. The current bioclimatic data (Biol 
to Bio19, for years 1950-2000) were downloaded from 
version 1.4 (30 arc-sec, or ~1 km) of the WorldClim 
website (http://worldclim.org, Hijmans et al. 2005). The 
varible definitions are given in Table 1. In Version 1.4, 
HadGEM2-ES (30 arc-sec, or ~1 km) based data for 19 
bioclimatic variables were available for the future (2050 
and 2070) climate projections based on the RCP 2.6, 
RCP 4.5, RCP 6.0, and RCP 8.5 scenarios. These data 
sets are provided on a global scale, and were optimized 
to the size of Turkey with the help of ArcGIS 10.2 and 
converted to ASCII format. 


Habitat suitability model. Due to the narrow study 
area, high correlations can exist between bioclimatic 


variables which may pose a problem during the analysis. 
To eliminate the problem of multi-collinearity, we 
applied Pearson Correlation Analysis (r? < 0.8) for the 
19 bioclimatic variables. If a pair of variables has a 
correlation coefficient of greater than 0.8, they were 
considered to both represent a similar phenomenon, 
so one of the pair of variables was excluded from the 
analysis. 

MaxEnt 3.4.1 software was used to estimate the 
climatic conditions limiting the current distribution of 
P. caralitanus and to compare the current situation with 
future climate scenarios (Phillips etal. 2017). The presence 
data of P. caralitanus in CSV format were entered into the 
“Samples” section of MaxEnt 3.4.1 software. The current 
bioclimatic data in ASCII format were entered into the 
“Environmental layers” section. Then, 2050 and 2070 
bioclimatic data files (RCP 2.6, RCP 4.5, RCP 6.0, and 
RCP 8.5) were entered in the “Projection layers directory/ 
file” section to compare the future scenarios. MaxEnt 
analysis was run using 90% of the records as training 
data to build the model and the remaining 10% for testing 
the model. Ten repetitions were made for each model so 
that different training and test data would be processed in 
the analysis. Any bioclimatic variables (Biol to Bio19) 
that did not contribute to the model obtained as a result 
of the analysis, were excluded in the next analysis. 
Replicated run type Crossvalidate was selected for the 
analysis along with the settings: maximum iterations 
500, convergence threshold 0.00001, log file maxent.log, 
and the default prevalence was set to 0.5. The analysis 
process continued until the best model was obtained. We 
used the “Area under the receiver operating characteristic 
(ROC) curves” (AUC) values to determine the best 


Table 1. Bioclimatic variables obtained from the WorldClim website (http://worldclim.org). 


Code Bioclimatic variables Unit 
Biol Annual Mean Temperature eC 
Bio2 Mean Diurnal Range (Mean of monthly (max temp - min temp)) aC 
Bio3 Isothermality ((Bio2/Bio7) * 100) unitless 
Bio4 Temperature Seasonality (standard deviation *100) CofV 
Bio5 Max Temperature of Warmest Month SC 
Bio6 Min Temperature of Coldest Month $C 
Bio7 Temperature Annual Range (Bio5-Bio6) °C 
Bio8 Mean Temperature of Wettest Quarter ne: 
Bio9 Mean Temperature of Driest Quarter “© 
Biol0 Mean Temperature of Warmest Ouarter 3G, 
Bioll Mean Temperature of Coldest Quarter = 
Biol2 Annual Precipitation mm 
Biol3 Precipitation of Wettest Month mm 
Biol4 Precipitation of Driest Month mm 
Biol5 Precipitation Seasonality (Coefficient of Variation) CofV 
Biol6 Precipitation of Wettest Quarter mm 
Biol7 Precipitation of Driest Quarter mm 
Biol8 Precipitation of Warmest Quarter mm 
Biol9 Precipitation of Coldest Quarter mm 
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Fig. 2. Variables with the highest contributions to the potential distribution of P. caralitanus according to MAXENT with the 
standard errors in blue. The Y-axis indicates the probability of presence (based on the Cloglog, or complementary log-log transform, 


values) and the X-axis shows the contribution of each variable. 


model performance. An AUC value close to 1 indicates 
that the model is in a perfect performance and a value 
of 0.7 or higher indicates that the model is descriptive, 
while a value of 0.5 indicates that the model does not 
provide useful information (Phillips et al. 2006). Among 
the models with the highest AUC values, the one with the 
lowest standard deviation was selected as the best model. 


Results 
The results of the analysis indicate that the AUC value 


of the training data is 0.965, while the test data AUC 
value is 0.964, and the standard deviation is 0.012. We 
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found that the climatic factors limiting the distribution of 
P. caralitanus were: Bio4 (percent contribution: 22.3%, 
Temperature Seasonality), Bio6 (16.4%, Min Temperature 
of Coldest Month), Bio7 (17.9%, Temperature Annual 
Range), Bio8 (19.9%, Mean Temperature of Wettest 
Quarter), Biol4 (15.5%, Precipitation of Driest Month), 
and Biol7 (8%, Precipitation of Driest Quarter). 
According to these results, P. caralitanus prefers habitats 
with precipitation in the driest months, and those where 
the average temperature of the coldest three months does 
not fall below -5 °C, the seasonal temperature difference 
is high, and the humid months are 0—5 °C (Fig. 2). 

The habitats for P. caralitanus were primarily in 
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Fig. 3. Current climatic habitat suitability map (A) and the eight RCP climatic change scenario maps for P. caralitanus based on 
RCP 2.6 (B-C), RCP 4.5 (D-E), RCP 6.0 (F—G), and RCP 8.5 (H-D for either 2050 or 2070 as indicated. 


Southwestern Anatolia, in the region between Central 
Anatolia and Western Anatolia, and around the western 
Taurus Mountains (Fig. 3A). A large part of the map 
reflects the distribution areas of P. caralitanus reported 
in the literature. Although its presence has not yet been 
reported from the new potential areas, these areas are 
climatically potential habitats. 

If the RCP 2.6 scenario occurs, there will be habitat 
losses in the Turkish Lake District until 2050, potential 
habitats in the Southwest will persist, and new potentially 
suitable bioclimatic habitats will emerge in the North 
Aegean region (Fig. 3B). By 2070, the situation in the 
Turkish Lake District will be slightly better than in 2050 
and new potential habitats will be formed in a small area 
in the inner Western Aegean region (Fig. 3C). In the case 
of the RCP 4.5 scenario, suitable habitats in the Turkish 
Lake District will almost disappear, new habitats will 
appear in the Northern Aegean regions, and potential 
habitats in the Southwest will persist (Fig. 3D—E). The 
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RCP 6.0 scenario map shows that a situation similar to 
the RCP 4.5 scenario in Fig. 3E will occur in 2050, but 
it reveals that the situation will be worse in 2070 (Fig. 
3F-G). According to the RCP 8.5 scenario, which is the 
most dramatic result, potential habitats in Southwestern 
Anatolia will be reduced but some will continue to exist, 
and the appropriate habitats of P. caralitanus around the 
Turkish Lake District will completely disappear. On the 
other hand, new potential habitats not seen in the other 
scenarios will emerge in the Northwest Aegean region 
(Fig. 3H-I). 


Discussion 


Recent studies have provided important information on 
the extinctions of species due to climate change. Studies 
have reported that there could be a 15—37% loss of species 
by 2050 according to various climate change scenarios 
(Thomas et al. 2004). Even in the most optimistic climate 
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change scenario, most of these losses will occur in 
endemic species (Dirnböck et al. 2011). In response to 
climate change, many species will present compensating 
mechanisms such as adaptation and migration; however, 
the success of these mechanisms against climate change 
will be related to the climate change velocity. In contrast 
to species which have large areas of distribution, 
endemic species distributed in a narrow region according 
to special climatic and topographic conditions will not be 
as successful, except in areas with a low climate change 
velocity. In areas with a high climate change velocity, 
endemism is either very low or does not exist. For this 
reason low climate change velocity areas will be essential 
refuges for many small-ranged species, such as endemic 
species (Sandel et al. 2011). 

The current habitat suitability map (Fig. 3A) shows 
the distribution of P. caralitanus mainly in the Turkish 
Lake District region and is compatible with species 
observation data. It also indicates that there are suitable 
habitats for this species in the north-western parts of the 
Turkish Lake District (Fig. 3A), while recent studies have 
found that the most north-western distribution area of the 
species is Acıgöl and Işıklı Lakes (Baskale et al. 2017). 
These lakes are located at the north-western border of 
the Turkish Lake District. Although there are locations 
with similar climatic conditions, the distribution of the 
species in the northern part is geographically isolated due 
to mountains. 

In the past two decades, some lakes in the Turkish Lake 
District have reduced water levels due to climatic changes 
(Göncü et al. 2017; Kantarcı 2008). Wetlands and lakes 
are the essential elements of the habitats of frog species 
such as P. caralitanus because they need water bodies to 
complete their life cycle. A significant reduction of snow 
in higher regions due to climate change, a reduction in 
the water storage basins and earlier melting of snow 
will have negative consequences for many wetlands in 
Turkey (Yilmaz et al. 2019). Potential distribution maps 
determined for P. caralitanus can help in prioritizing the 
planning of future wetland use management around the 
current populations, the discovery of new populations or 
the identification of top priority areas. 

According to RCP scenarios that were created here for 
P. caralitanus, the climate change velocity will be high 
for the next 30 to 50 years (Fig. 3B—I), and habitat loss 
and fragmentation will be quite high in the distribution 
area of the species. For the survival of endemic species, 
we found RCP 2.6 to be the most optimistic scenario 
among the climate scenarios. However, even in this 
scenario, if P caralitanus cannot migrate to the north- 
western Turkish Lake District, the distribution of the 
species will shrink in the next 30 years, and a very 
narrow distribution area will host the species in 2070. 
The IPCC, which is currently in its 6" assessment phase, 
reported “Global warming of 1.5 °C” towards the end 
of 2018, and emphasized that a temperature rise by 1.5 
°C was better than 2 °C (IPCC 2018). According to our 
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results, if the RCP 2.6 scenario is exceeded, most small- 
scale wetlands will begin to dry in the next 30 years, and 
irreversible changes may occur in the following years 
relative to the current distribution map of P. caralitanus 
(Fig. 3D-I). Although new locations with similar climate 
zones occur in the north and north-west of the Turkish 
Lake District, it would be very difficult for the frogs to 
migrate from the original habitats to the new climatically 
suitable habitats within the next 30 and 50 years because 
of geographic isolation. Moreover, if the RCP 8.5 climate 
scenario occurs, there will be no climatically suitable 
habitats for P. caralitanus in the Turkish Lake District, 
which we can categorize as a disaster. 

We have shown that the MaxEnt model is a useful tool 
for creating a predictive distribution map of P. caralitanus 
in the Turkish Lake District, and it has allowed us 
to determine putative environmental constraints and 
successfully predict the species potential fate in the 
coming years. According to our results, P. caralitanus 
will move farther towards extinction in future, in the 
face of such problems as an increase in temperature, 
decrease in precipitation, loss of habitat, and reduction of 
water bodies. In order to ensure the continuation of this 
endemic species, the protection of habitat and the impact 
of climate change on the species should be investigated 
in detail with more environmentalist approaches. 
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Appendix 1. The list of localities of Pelophylax caralitanus in this survey. 


Site No 


O AN DN FW NY — 


peak peek 
a il) 


45 


Climate change and Pelophylax caralitanus 


Locality 


Acıgöl site 1/Denizli 
Acıgöl site 2/Denizli 
Acıgöl site 3/Denizli 
Acıgöl site 4/Denizli 
Acıgöl site 5/Denizli 
Ağlasun site 1/Burdur 
Ağlasun site 2/Burdur 
Ağlasun site 3/Burdur 
Akgöl site 1/Burdur 
Akgöl site 2/Burdur 
Akgöl site 3/Burdur 
Aksu river/Antalya 
Akşehir lake site 1/Konya 
Akşehir lake site 2/Konya 
Akşehir lake site 3/Konya 
Belceğiz/Isparta 
Beyşehir site 1/Konya 
Beyşehir site 2 /Konya 
Beyşehir Site 3/Konya 
Beyşehir site 4/Konya 
Beyşehir site 5/Konya 
Bucak/Denizli 
Burdur lake site 1/Burdur 
Burdur lake site 2/Burdur 
Burdur lake site 3/Burdur 
Derebucak site 1/Konya 
Derebucak site 2/Konya 
Eber lake site 1/Afyon 
Eber lake site 2/Afyon 
Eber lake site 3/Afyon 
Eber lake site 4/Afyon 
Eğirdir lake site 1/Isparta 
Eğirdir lake site 2/Isparta 
Eğirdir lake site 3/Isparta 
Eğirdir lake site 4/Isparta 
Eğirdir lake site 5/Isparta 
Eğirdir lake site 6/Isparta 
Fele/Isparta 
Gedikli/Isparta 
Gencek/Konya 
Gökgöl/Denizli 
Gökhöyük/Konya 
Gölcük site 1/Isparta 
Gölcük site 2/Isparta 
Gölcük site 3/Isparta 
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GPS Coordinates 
Latitute Longitute 
37°46.375°N 29°50.680’°E 
37°49.794 N 29°56.879 E 
37°46.348°N 29°50.58 P E 
37°47.003’N 2951 709'E 
37°50.942’N 29°59.452’E 
37°36.065’N 30932.018'E 
37935.967'N 30932 420'E 
37936.195'N 30932 406'E 
37°41.323N 29°44.759°E 
37°40.268 N 29°44 206'E 
37°39.887 N 29945 293'E 
37914.896'N 30947. 43VE 
38926.834'N 31°21.994’E 
38°32.596’N 31°22.83F E 
38°31.31# N 31°19.875 E 
37°57.680°N 31°19.188 E 
37958.268N 31925.98/'E 
31233332 N 31°30.742’E 
37°47.691°N 31°35.816 E 
37°39.479 N 31°41.819 E 
37°44.304°N 31°41.056’E 
38915.035'N 29°51.89P E 
37°41.411’N 30903.452'E 
37938.995'N 30903.628'E 
3738225 N 30°05.685 E 
37°21.654°N 31932.842'E 
37925.667'N 31°30.63 P E 
38°36.832’N 31S: 567E 
38°39.140 N 31°14.51S°E 
38°42.375°N 31°11.50 PE 
38937.808N 31904.053'E 
38°05.178 N 30°56.857 E 
38°16.532’N 30951.558VE 
38°16.723 N 30°52.190 E 
38°01.23 P N 30948.872'E 
38°08.23 N 30°52.789 E 
3700550 N SÖZE 
38°02.435’N 3177455 E 
37953.334'N 51205921 
326572 N 31°32.438°E 
38°12.207 N 30°02.938 E 
37924.056'N 31°56.421’E 
37°43.707 N 30°29.510 E 
37°44.469 N 30°29.030 E 
37°44.469 N 30°29.030°E 


84 
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Altitude (m) 


845 
84] 
840 
839 
840 

1,024 

1,023 

1,024 
994 
992 
994 
103 
954 
954 
954 

7 

1,132 

1,126 

1,136 

1,128 

1,126 
822 
910 
859 
865 

1,261 

1,281 
971 
970 
970 
970 
920 
921 
921 
921 
921 
921 

1,234 

1,127 

1,354 
821 

1,021 

1,387 

1,346 

1,347 
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Appendix 1 (continued). The list of localities of Pelophylax caralitanus in this survey. 


Site No 


46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
17 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 


Locality 


Gölhisar channel/Burdur 
Gölhisar dam/Burdur 
Gölhisar Lake/Burdur 

Gölyaka/Konya 
Gölyüzü/Konya 
Huğla/Konya 

Işıklı lake site 1/Denizli 

Işıklı lake site 2/Denizli 

Işıklı lake site 3/Denizli 

Işıklı lake site 4/Denizli 

Işıklı lake site 5/Denizli 

Işıklı lake site 6/Denizli 

Işıklı/Denizli 
Karaköy site 1/Burdur 
Karaköy site 2/Burdur 
Karamık site 1/Afyon 
Karamık site 2/Afyon 
Karamık site 3/Afyon 

Karaöz/Antalya 

Kavakköy/Denizli 

Kırkgöz site 1/Antalya 
Kırkgöz site 2/Antalya 
Kırkgöz site 3/Antalya 

Kovada/Isparta 
Kumluca site 1/Burdur 
Kumluca site 2/Burdur 

Salda lake site 1/Burdur 

Salda lake site 2/Burdur 

Salda lake site 3/Burdur 

Salda lake site 4/Burdur 

Salda lake site 5/Burdur 

Sarayköy/Konya 
Sarıkaya/Isparta 
Seydişehir site 1/Konya 
Seydişehir site 2/Konya 
Seydişehir site 3/Konya 
Sugözü/Isparta 
Suğla lake site 1/Konya 
Suğla lake site 2/Konya 
Suğla lake site 3/Konya 
Taşağıl/Konya 
Tepearası/Konya 
Yalıhüyük/Konya 
Yeşildağ/Konya 
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GPS Coordinates 
Latitute Longitute 
37°02.734 N 29°27.960 E 
37°01.385’N 29°26.687°E 
37906.71107N 20X55. 370'E 
37°41.825°N 31926.898'E 
37221-225 N 31°55.42F E 
3T 28112 N 31°34.050 E 
38°11.01# N 29°46.840 E 
38°15.332 N 29°51 .81 PE 
38°11.176 N 29°59. 55P E 
38914.403'N 29957245 
38916.443'N 29954 415'E 
38°12.368 N 30°01.774E 
38°19.25S’°N 299513751E 
37°28.333 N ETTE 
37°28.56 N 29932 385 E 
38°27. 915°N 30°52.449°E 
38°23.005’N 30°45.39VE 
3820175 N 30°53.148°E 
Sie O67 N 30948.640'E 
38909 .254'N 29°38.27 P E 
37°06.600’N 30934 828'E 
37°05.669 N 30535252E 
37°05.526 N 30°34.997 E 
37230 7534N 30953 886'E 
37°38.270 N 30°02.836 E 
37°38.518 N 30°03.275F E 
37°31.39P N 29°38.159 E 
37°31.142’N 29°37.493°E 
37441725. N 29°39. 41VE 
37318185 N 29937. 40'E 
37°31.499 N 29936 012 E 
37919.247'N SOME 
37°56.404°N 31°18.464E 
37°28.588 N 31°48.769 E 
37°28.519 N 31949 462'E 
37°27.094 N 31°51.63P E 
38°14.227 N 31°17.23P E 
37°21.50 N 31°59.417 E 
37918.965'N 32°06.145’E 
37221232N 31255510 E 
37°22:702:ŅN 31953321E 
37°28.233 N SISKS at 
37919.256'N 32906.28/'E 
37°35.087 N 31°30. 1SVE 
85 
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Altitude (m) 


1,022 
1,068 
947 
1,126 
1,099 
1,370 
816 
819 
821 
821 
821 
821 
829 
996 
991 
1,004 
1,004 
1,004 
55 
811 
304 
307 
302 
906 
860 
860 
1,158 
1,161 
1,145 
1,160 
1,168 
1,102 
1,127 
1,135 
1,131 
1,110 
1,226 
1,222 
1,096 
1,098 
1,101 
1,387 
1,101 
1,125 


